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Abstract 

The vortex state thermal and transport properties of the high Tc copper ox- 
ides can be understood in a d-wave gap model and are dominated by the ex- 
tended quasiparticle states that exist along the nodal directions in momentum 
space. The Doppler shift on these states due to the circulating supercurrents 
around the vortex core, introduces new van Hove ridges into the energy de- 
pendent local density of states (LDOS) as a function of distance in the region 
between cores. We emphasize the topology of these ridges and the effect on 
them of local static impurity scattering in Born and unitary limit. We treat 

possible orthorhombicity. Effective local scattering rates are also obtained. 
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I. INTRODUCTION 

The theoretical predictions of Volovik0i and Kopnin and Volovikifl that the low temper- 
ature specific heat in a (i-wave superconductor would vary as the square root of the applied 
magnetic field {^/H) was verified in several experiment s-iHi To understand the experimental 
data it is assumed that the extended quasiparticle states outside the vortex cores can be 
treated semiclassically taking account only of the Doppler shift due to the circulating super- 
currents. To obtain the specific heat, which is a bulk property of the system, it is necessary 
to carry out a spacial average over the vortex unit cell of the local density of states (LDOS) 
which varies as a function of distance from the vortex core center. The low temperature 
limit of the specific heat depends on the zero frequency limit of this averaged density of 
quasiparticle states. This quantity is not only affected by the magnetic field H but can 
also depend on impurity content ,00 particularly in the resonant scattering limit. In this 
limit the theory of a rf-wave superconductor without magnetic field predicts a finite, con- 
stant value for its zero frequency limit which is proportional to the self consistent effective 
scattering rate and this leads to the well known universal limit for the transport .SHll 

When the magnetic field is oriented in the copper oxide plane rather than perpendicular 
to it new anisotropy effects in the specific heat are predicted.0'llHli In addition transport 
properties are also affected by the Doppler shift of a vortex statcllllli The thermal conduc- 
tivity has been particularly widely studied. 

The semiclassical Doppler shift of the electron due to the circulating supercurrents out- 
side the vortex core should also affect the frequency dependence of the LDOS which can 
be measured in scanning tunneling microscope (STM) experiments. In view of the above 
successes in our understanding of the thermal and transport properties of the copper ox- 
ides it is of interest to understand the signature of the Doppler shift on the LDOS and 
the effect of impurities on it. Very recently Franz and TesanovicH considered the LDOS 
problem and compared an average over vortex winding angles of the semi classical LDOS 
with equivalent results obtained from complete solutions of the BdG equations for a single 
vortex. Their main conclusion was that the semiclassical approach does indeed provide a 
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good approximation for this quantity. 

For no magnetic field {H = 0) there is a single van Hove singularity in the density of states 
at energy u equal to the gap amplitude which does not shift with position r (homogeneous 
case). For finite H new van Hove ridges are predicted. The actual number of such ridges 
and their topology depends on the winding angel (3 and simple formulas can be obtained 
from which their contour in (cj, r)-space can be traced. The effect of impurities on the 
topology of these ridges is considered as is the effect of orthorhombicity. Impurities are 
treated in both Born and unitary limit, and the influence of the circulating supercurrents on 
the local effective impurity scattering is also considered. This topic was recently discussed 
by Barash and Svidzinskii.0 Orthorhombicity, present in YBa2Cu3 06.95 (YBCO) because 
of the existence of chains, is treated within a simple effective mass model for the in-plane 
electronic dispersion curves and a possible subdominant s-wave component is added to the 
dominant c/-wave gap. 

The tetragonal case is the subject of section II. Impurities are included. Generalization 
to the orthorhombic case is given in section HI. Section IV deals with effective local impurity 
scattering rates. A short conclusion is found in section V. 



II. TETRAGONAL CASE 

In a semiclassical approximation the local Green's function at position r is given in terms 
of the Doppler shift Vi7'(k)qs by 

G{k, Un] r) = — 2 , 2 I A2 (l) 

where ti^2,3 are the 2x2 Pauli matrices, iuon the Matsubara frequencies i{2n + l)'7rT, n = 
0, ±1, ±2, . . . and T is the temperature. The electronic dispersion in momentum space k is 
and the gap in the two-dimensional copper oxide Brillouin zone is A^. The electronic 
Fermi velocity is Vj7(k) and is the momentum associated with the superfiuid flow about 
the vortex. In the magnetic field region Hd < |H| ^ Hc2 with ilci,2 the lower and upper 
critical field respectively, and H the external magnetic field taken to be perpendicular to the 
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Cu02-planes, the intervortex distance R = ^y^^y where a is a geometrical factor of order 
one which is associated with the vortex arrangement, and 0o is the fundamental quantum 
of flux. In what follows, the position r will be measured relative to the center of the vortex 
core. We use polar coordinates {r, (3) with /3 the vortex winding angle. The dimensionless 
variable p = r/R takes on the value 1 at the boundary of the vortex unit cell. It is expected 
that R^ C,o (^0 is the coherence length) and we will be interested only in values of r > ,^o- 

Assuming a circular velocity fleld for the supercurrents around a single vortex, = 
with /3 a unit vector along the current, the Doppler shift 

vp(k)q, = ^sin(0-/3), (2) 

where Eh is the magnetic energy scale that enters our problem and is equal to ^vp^J^- = 
uAq where the last identity measures Eh in units of the gap amplitude Aq. In Eq. (|^) </> 
is the polar angle associated with momentum k on the Fermi surface and the d-wave gap 
Ak = Aocos2(/) = in the same model. 

The local quasiparticle density of states at relative position p and energy u, Nl{p,u!), 
follows directly from the analytic continuation of G deflned in Eq. ([l|) to real frequencies 
iujn —>■ uj + iO^ with 

^^^^'^^= /^g?e| , l"-"Hk)q.| 1 







- v^(k)q,]^ - A2(0) 



where Nq is the normal state density of states. Note that the Doppler shift of Eq. @ which 
enters depends on the magnetic fleld |H| = if through Eh ~ VH and on the position 
of the STM tip away from the vortex core given by p and (3. We will be interested only 
in the region where r is outside the vortex core of size ^o- This is mandated by the fact 
that our approach does not treat in detail the core interior but instead includes only the 
supercurrents outside the core through the superfluid velocity fleld = which itself is 
an approximation. Nevertheless, it is the aim of this work to understand how the presence 
of Vs modifles the LDOS in the region between the vortex cores. 

The integral in Eq. (|^) is easily evaluated. We present a flrst set of results which in- 
clude some impurity scattering in Born approximation. This smooths out the logarithmic 



singularity at the gap amplitude predicted for pure d-wave with H = 0. To include static 
impurity scattering it is simply necessary to exchange the frequency u and the gap Arf(0) in 
Eq. (^) by their renormalized values denoted by a tilde, i.e.: oj and A(0). The renormalized 
frequency 

u{p, l3,uj) = uj + iTTt^il{p, /3,uj) (4) 
in Born approximation with t+ the impurity scattering rate in the normal state, and 

Uj{p,P,Uj} = UJ + tTTT — r (5) 

in the unitary limit with F"*" replacing t"*" of Eq. (p. The renormalized gap A{(f)) = 
Ao cos 2(j) + iAs{p, f3, u) in a (d + zs)-wave symmetry with 

Asip,/3,uj) = tTTt+D{p,/3,iu) (6) 

for Born scattering and 

A, p, f3, u) = «vrr+— — — — . . — r (7) 
n^{p,f3,uj) + D^{p,f3,u;) 

in the unitary limit. Here, the two functions Q{p,P,u!) and D{p,P,u!) are 

\ ^Hp, P, ^) - vr(p)q,l^ - A^(p, P, ^) - A^(0) 

and 

£'(p,/?,^) = ( ^ ) , (9) 

\ ^^(P, /?, a;) - vp(p)q,]^ - A2(p, /?, - A2(0) / ^ 

where (■ ■ ■)0 indicates an angular average over the angles around the circular Fermi surface 
in the two-dimensional copper oxide Brillouin zone. 

In the top frame of Fig. |l] we show our results for a d-wave superconductor (tetragonal 
symmetry) with a gap amplitude Aq = 24meV. We take the magnetic field H || c-axis, i.e. 
H is perpendicular to the copper oxygen planes and the vortex winding angle /? = 0°. The 
vertical axis is labeled by Nl{p,uj) which is given by 
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iV^(p,a;) = 3?e{n(p,Aa.)}^_,, (10) 

and features the LDOS as a function of p and uo. Estimates of the magnetic energy scale 
Eh in YBCO have been considered by Vekhter et al.0 This quantity is not very well known 
but best estimates give between 



Eh = 20 - 30[K]^i7[Tesla]. 

Therefore the magnetic energy scale has been set at 10% of the gap amplitude for illustrative 
purposes, i.e. z/ = 0.1 in our notation and the static impurity content in Born approximation 
is = 0.1 meV, enough to smear out the van Hove singularity but not enough to make 
them disappear. (This value of corresponds to a scattering rate of about 6 K. Recent 
experiments by Hosseini et al.0 seem to indicate an even smaller scattering rate.) Three 
prominent sets of ridges are seen for p in the range 0.1 (at the vortex core, i.e. at r = ^o) and 
1.0 (at the vortex cell boundary). The case with no magnetic field is shown as the bottom 
frame of Fig. |l| for comparison. Only one prominent van Hove ridge is seen in this second 
case and it is centered at u = Aq = 24meV and, of course, there is no p dependence since 
there is no vortex core and thus no spatial inhomogeneities in the problem. A main ridge 
around u = 24 meV is still seen in the top frame although as p is decreased and the vortex 
core is approached at p = 0.1 it has shifted in energy significantly. The other two ridges 
are due entirely to the Doppler shift provided by the supercurrents and can be understood 
as follows. In the pure case (no impurity scattering) the ridges in Nl{p,uj) correspond to 
extrema in the zeros of the denominator of Eq. (|^), i.e. they correspond to the extrema of 
the equation 



uj = v^(k)q, ± A(0) = Ao 



— sin(0 — (3) ± cos 20 
P 



(11) 



i.e. to solutions of equation 



= = -cos(0-/?) T2sin20, (12) 
d(f) p 

which have been written for a general value of vortex winding angle (3. There are three 
solutions of this equation for the case considered in the top frame of Fig. |I] (/3 = 0). The 
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critical values of frequency ujd with i = 1,2,3 are shown as a function of distance from the 
vortex core p in the top frame of Fig. ^. For u > the solutions occur at the critical angles 
0ci,2 = ±vr/2 and (pcs = sin~"'^(z//4p) respectively and the corresponding values of tOci are 

a;a,2 = Ao(l±^^/p), (13) 
which are valid without restriction; in the limit u/p ^ we find 

2" 



2 l2p, 



(14) 



Numerical solutions can always be obtained and are needed except for a few simple cases. 
For 1^ = 0, no magnetic field, all Ud = Aq as expected. At the vortex cell boundary p = 1, 
i.e. r = R, the three critical frequencies are in units of Aq 1.1, 0.9, and 1.0013 respectively, 
and at the vortex core which corresponds to p = 0.1 for u = 0.1 these values are 2, 0, and 
1.125. These expectations are bourn out in the top frame of Fig. |I|. We note from Eq. (13 



that the spacing of the two additional structures away from u = Aq goes in the antinodal 
direction (/3 = 0°) like 



Ucl,2 (X VH, 

i.e. is linear in u; the shift of the main van Hove ridge from Aq due to the Doppler shift 
goes instead like 

Uc3 oc H, 

i.e. is linear in the magnetic field in the limit z//p — 0. Other directions, (3 ^ 0°, will 
develop other if-dependencies of the spacing of the van Hove ridges which will have to be 
evaluated numerically in most cases. 

Also note the — limit of the LDOS. It goes to zero linearly in the H = case but, in 
contrast, it takes on a finite value when vortices are present. Its value depends strongly on 
position p in the vortex unit cell and increases strongly with decreasing p. It is the average 
of the LDOS at = over the vortex unit cell, i.e. the average over p and winding angle (3 
that determines the low temperature bulk specific heat and this varies with the value of H. 



Franz and Tesanoviccj discuss in detail how this region of small uj can be strongly modified 
by an in-plane angular dependence to the tunneling matrix element that comes into the 
relationship between STM measurements and the LDOS itself. These are only the same 
when the tunneling element is constant. 

In the bottom frame of Fig. Q we show again the same data as shown in the top frame 
of Fig. [T| but in a different representation which allows to see the ridges a little more clearly. 
What is presented here is Nl{p^uj) vs. uj for fixed values of p, namely p = 0.2 (solid line), 
p = 0.3 (dashed line), p = 0.6 (dotted line), and p = 1.0 (dash-dotted line). The vortex core 
radius equal to the coherence length corresponds to p = 0.1 and the inter vortex distance 
is p = 1.0. Finally, a thin solid line indicates Nl{uj) for no magnetic field, H = u = 0. 

Furthermore, the data shows clearly the main van Hove singularity around, but not 
precisely a.t uj ^ Aq, with some variation with p as detailed in the top frame of Fig. |^ 
(solid line, ujcs). This peak exists even when H = (light solid line) when there are no 
supercurrents and it is the usual rf-wave gap peak. Note that here we have included some 
small amount of impurity scattering in Born approximation, namely = 0.1 meV. The 
other secondary peaks above and below u ~ 24meV are due to the supercurrents and are 
the semiclassical signature of the resulting Doppler shifts. (Similar ridges with somewhat 
different geometry appear in earlier work which makes use of the semiclassical Eilenberger 
approach.&H This approach should be better inside the core, a region we do not treat 
properly. Other related work in the s-wave case is the very recent paper by Eschrig et al.il) 
As the STM tip is moved more towards the vortex core, the secondary peak at low energy 
becomes fairly prominent and also the u = limit becomes finite as noted before. The 
topology of the surface Nl{p,uj) vs. p and u is affected by several external variables such 
as the magnitude of the magnetic field, the vortex winding angle, the type of impurities 
involved, i.e. Born or unitary scattering, and its strength as we now detail. 

The ridge structure in the LDOS changes considerably as the position of the STM tip 
is changed. This is illustrated in the top frame of Fig. ^ where we show our result for a 
winding angle (3 = 45°, i.e. in the nodal direction. There are striking differences with the 
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top frame of Fig. |l[ Only two ridges are now seen in contrast to three in the (3 = 0° case 
(the antinodal direction), and the main gap peak around u = 24meV at the vortex cell 
boundary p = 1 {at r = R) moves downward in energy as p is reduced toward the vortex 
core. The a; = value of Nl{p, uj) near p = 0.1, i.e. r = C,o, is also much smaller than in the 
(3 = 0° case. 

In the bottom frame if Fig. ^ we show our solutions for the critical values of uj {ujci, i = 
1,2) which locate the ridges in the top frame. The critical angles ipdii = 1,2 at which the 
extrema responsible for the ridges occur are also presented as functions of the distance from 
the vortex center p in this frame. It is obvious that geometry strongly affects the topology 
of the Nl{p,lj) surface. 

Impurity scattering can also affect things in an important way. In Fig. ^ we show results 
for the LDOS Nl{p,uj) as a function of the distance to the vortex core p and energy uj for 
a magnetic energy z/ = 0.1, H || c-axis, and a winding angle (3 = 0° (antinodal direction). 
This figure is to be compared directly with the top frame of Fig. ^ as the only difference 
between the two figures is that in Fig. ^ the unitary impurity limit is used, Eqs. (^) and 
(0), with = 0.1 meV instead of the Born limit, Eqs. (^) and (^), with = 0.1 meV. The 
structure of the ridges is now more pronounced and the zero frequency limit of the LDOS 
near the vortex cell boundary is now larger. 



III. ORTHORHOMBICITY 

YBCO is orthorhombic because of the existence of chains along the 6-axis in one of the 
three copper oxide two-dimensional planes per unit cell. This implies that the gap is in 
principle a mixture of s- and (i-wave symmetry which do belong to the same irreducible 
representation of the crystal lattice group, and therefore has the form 

A(0) = Ao(cos20 + s) (15) 

on the Fermi surface. Here s is the subdominant s-wave component. A simple way to include 
band anisotropy is to introduce anisotropic effective masseJiil@ with < nia in an infinite 
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band model. This gives an ellipsoidal Fermi surface which can always be mapped onto a 
circular one by scaling of momentum variables. At the same time the gap is transformed 



A(0) = Ao s + 



a + cos 2(j) 



(16) 



1 + a cos 20 ) 

There are no additional changes and the calculations proceed as before with Eqs. (pA|) and 



(|T^). For the direction /3 = 0°, they are replaced by 

= Ao 

and 



p 1 + a , / a + cos 20 

sm (p ± I z — + s 



p VI - a 



1 + a cos 20 



(17) 



duo 







z/ 1 + a , 2(1 - a^) sin! 

: COS =F 



aip p \J\ — a ^ ' (1 + a COS 20)2' 

with a the effective mass anisotropy parameter defined as0S a = {nib — '"^a)/(^a + rrib)- 
Again, for a general winding angle numerical solutions are needed although analytic ones 
are still possible in the limit u/p —>■ in some cases. If we restrict ourselves to solutions of 
Eq. (|16D which have at least some positive values for the critical frequency Uc in the range 
z/ < p < 1 then we get for j3 = 0°, and the critical angles 0c = ±7r/2 the following three 
solutions of Eq. (|l] 



a; 



i±^±(l-s), 

+ (1 - ^) 



p \/l—a 
V l+g 



(19) 



p \J\—a 



without restrictions on the value of z//p, and another solution 



1 + s 



a 



Ao ■ ■ " ■ 2 V2p/ (1 



a] 



(20) 



which is valid only for v / p and corresponds in that case to 0c near zero radians. 

First we note the well known limit of z/ = (no magnetic field) in which case the van 
Hove singularity in the density of states is split by the existence of the subdominant s-wave 
component of the gap and there are two singularities in the density of states positioned at 
ujc = Ao(l ± s). This holds whatever the value of the effective mass anisotropy a might be. 
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This mass anisotropy does not lead to a corresponding splitting of the van Hove singularities 
in the density of states. In a previous paper Schiirrer et al.@ have considered the specific 
case of optimally doped YBCO and have suggested that reasonable model parameters are 
a = 0.4 to account for a factor of a little more than 2 in the effective conductivity between 
b- and a-direction. Further consideration of finite temperature penetration depth data@ 
suggests s = —0.25. This gives agreement between the model and the measured penetration 
depth data at low but finite temperatures T. 

In the top frame of Fig. ^ we show results for Nl{p, uj) for a case with a = 0.4, s = —0.25, 
and f3 = 0°, i.e. in the antinodal direction. The top frame is a three dimensional plot of 
Nl{p,uj) vs. p and u and is to be compared with the top frame of Fig. [l| (tetragonal case). 
We note striking differences in the ridge structure. Now, the two ridges at p = 1 and about 
UJ = Ao(l±s) come together as the vortex core is approached and they eventually cross each 
other before diverging in the opposite direction. The contours for the critical frequencies, 
ujc, are shown in the bottom frame of Fig. ^. The two ridges that cross are the solutions of 
Eq. (|I^) for 0c = — 7r/2 and of Eq. (^^. The solution of Eq. ( [T^ ) for 0^ = 7r/2 and the plus 
sign gives the higher energy ridge shown in the top frame of Fig. |^ (and dash-dotted line in 
the bottom frame of this figure). 

As in the tetragonal case (Fig. ^ the ridge topology in the LDOS changes considerably 
as the position of the STM tip is changed. We illustrate this in the top frame of Fig. |^ 
where we present results for a winding angle (3 = 45°. Now all four ridges are seen in the 
entire range < p < 1 and two cross overs can be observed. The contours of the critical 
frequencies Uc, are shown in the bottom frame of Fig. |^. A comparison of the top frame 
of this figure with the top frame of Fig. ^ (tetragonal case, /3 = 45°) makes the splitting 
of the van Hove singularity as a result of the subdominant s-wave component of the gap 
particularly transparent. 
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IV. LOCAL IMPURITY SCATTERING RATE 



Finally, it is of interest to discuss briefly the local impurity scattering rates which are 
modifled by the presence of supercurrents. Similar work has already appeared in a paper by 
Barash and Svidzinskii.0 In Born approximation '^mu{p, P,uj) given by Eq. @ is given by 

Qmu{p,P,uj) =7rt^^en{p,P,uj), (21) 

where ^eQ{p, P,uj) is just the LDOS of Eq. (^ modifled by the impurity scattering and 
plotted already in the top frame of Fig. |I| for the antinodal direction and in Fig. ^ for 
the nodal direction. It is clear then that the Doppler shift accounting for the supercurrents 
signiflcantly affects the local scattering rates in the same way as they modify the LDOS. The 
changes due to the supercurrents are very different when unitary scattering is considered. 
In this case, we flnd from Eq. (1), 

SmcD(p,/3,u;) = 7rr+a^?e/ ^^^f'^'^o! n ^ l• (22) 

For a flxed winding angle /3 deflning 7(p, u) = Qmu}{p, (3, u) we have for u = 0.1, H || c-axis, 
and P = 0° the results shown in the top frame of Fig. |^ where we give a three dimensional 
plot of 7(p, uj) as a function of distance p from the vortex core and of energy uj for a system 
of tetragonal symmetry. This quantity is strikingly different from its Born limit counterpart 
which just corresponds to the Nl{p,uj) as plotted in the top frame of Fig. |I] scaled by a 
constant factor of vrt"^. 

The H = case in the unitary limit is quite similar to the H ^ case at p = 1. The only 
difference is that for H = there is only one shallow dip at u = Aq = 24meV which does 
not change position as one moves from p = 1 to p = u. It is quite important to notice that 
in the limit uj —>■ '~f{p,uj) decreases by almost one order of magnitude as one approaches 
the vortex core. This means that the supercurrents surrounding the vortex core will also 
influence substantially the value of the universal limit for the transport. 

The equivalent results for a system with orthorhombic symmetry with an effective mass 
anisotropy parameter a = 0.4 and an s-wave component s = —0.25 to the gap are shown in 
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the bottom frame of Fig. |^ to contrast the results for the system with tetragonal symmetry. 
The variation of '~f{p,uj) in the limit — is less dramatic in this particular case. 

It is obvious from these results that the supercurrents modify significantly the local 
effective impurity scattering which depends strongly on p. This has already been emphasized 
by Barash and SvidzinskiiH in another context. 

V. CONCLUSION 

In this work we have described the effect of supercurrents around the vortex core on 
the LDOS and on the local impurity scattering. We have included in the work only the 
effect of the Doppler shift in a semiclassical approximation on the extended quasiparticles. 
The work was motivated by the considerable success such an approach has had recently in 
describing the effect of vortex cores on the thermodynamics and some transport properties of 
(i-wave superconductors for magnetic fields in the range Hd < H <^ Hc2, between the lower 
and upper critical fields. The calculations predict definite modifications of the expected 
van Hove singularity at the gap amplitude Aq in a (i-wave superconductor. There are 
additional ridges which depend on the geometry of the STM arrangement, i.e. its distance 
from the vortex core and direction of the STM with respect to nodal or antinodal direction. 
They also depend on orthorhombicity and are modified by impurity scattering. So far we 
are not aware of experiments which confirm these predictions. We hope our work will 
stimulate more experiments. Of course more sophisticated approaches to the structure of 
a vortex core in a rf-wave superconductor are possible. We have already mentioned 
the work of Franz and Tesanovic.ll They compare results for the LDOS obtained from a 
semiclassical approach with results from solutions for a single vortex within a BdG self 
consistent approach. For the angular averages over the vortex winding angle they consider, 
they do find good agreement between the two sets of results. An important additional 
aspect of these author's work is that they investigate the possibility that the tunneling matrix 
element may have angular dependence and that consequently STM does not directly measure 
the LDOS. In the particular model considered, the perpendicular tunneling is modulated 
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by a cos^ 20 in-plane dependence which strongly affects the low energy part of the STM 
spectrum because there can be no tunneling right on the diagonal of the Brillouin zone. 
This modulation, however, will have less effect away from the a; ~ region and is not 
important for the ridges. Here our own emphasis has been on the topology of the secondary 
van Hove ridges introduced in the LDOS by supercurrents. We have derived simple formulas 
for the position of these ridges in configuration space, for their dependence on H which varies 
with position about the vortex, and have examined the effect of impurities both in Born 
and in unitary scattering hmit on the topology of the ridges. We have also considered the 
modifications that are introduced when orthorhombicity is introduced within an anisotropic 
effective mass model. 
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FIGURES 

FIG. 1. The top frame is the local density of states N^^p, to) as a function of energy iv and the 
distance from the vortex core p for a magnetic energy Eh = z^Aq with u = 0.1 and a vortex winding 
angle P = 0° along the antinode in the d-wave gap. A small amount of elastic impurity scattering 
in Born approximation has been included with = 0.1 meV. The main ridge around 24meV is the 
usual van Hove singularity at the gap amplitude Aq modified only slightly by the magnetic field 
H. The other two ridges are due to the Doppler shift of the semiclassical approximation resultant 
from the existence of the supercurrents about the vortex core. The bottom frame is the density of 
states when H = and is for comparison. 

FIG. 2. The top frame gives the critical values of frequencies coci, i = 1, 2, 3 defining the three 
ridges in the top frame of Fig. |l]. The bottom frame represents the same data as in the top frame of 
Fig. H but gives Ni{p,uj) vs. uj for several values of />, namely p = 0.2 (solid line), p = 0.3 (dashed 
line), p = 0.6 (dotted line) and p = 1 (dash-dotted line). Here, /o = 0.1 corresponds to the vortex 
core radius equal to the coherence length and p = 1.0 the vortex cell of Radius R. Included is 
also, for comparison, N{uj) for = 0, i.e. no magnetic field, as a thin solid line. 

FIG. 3. The top frame shows a three dimensional plot of the LDOS Nl{p,uj) as a function of 
distance from the vortex core p in the vortex unit cell, and of frequency uj. Here the vortex winding 
angle is /3 = 45°, i.e. the STM tip is moved along the nodal direction, and the magnetic energy is 
u = 0.1. The bottom frame gives the critical frequency uJci,i = 1,2 and the corresponding critical 
angles (j)ci that locate the ridges in the top frame as a function of the distance from the vortex core 
p, for u = 0.1, H II c-axis, /? = 45°. 

FIG. 4. A three dimensional plot of the LDOS as a function of the distance p from the vortex 
core in the vortex unit cell, and of frequency oj. Here the magnetic energy is = 0.1, H || c-axis, 
and the vortex winding angle /3 = 0° (antinodal direction). Unitary impurity scattering is included 
with r+ = 0.1 meV as compared to the top frame of Fig. |l| where the Born limit is used instead. 
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FIG. 5. Top frame: a three dimensional plot of the LDOS Nl{p,uj) as a function of distance 
from the vortex core p in the vortex unit cell, and of frequency co. Here the vortex winding angle 
f3 = (antinodal direction) and H || c-axis, i.e. perpendicular to the copper oxide planes. The mass 
anisotropy a = 0.4 and an s-wave component s = —0.25 are included to model the orthorhombicity 
of YBCO. Born impurity scattering with t+ = 0.1 meV is also included. The bottom frame gives 
the critical frequencies Uc defining the van Hove ridges in the top frame. 

FIG. 6. The top frame shows a three dimensional plot of the LDOS Nl{p,u) as a function of 
distance from the vortex core p in the vortex unit cell, and of frequency uj. Here the vortex winding 
angle j3 = 45° and H || c-axis. The effective mass anisotropy a = 0.4 and an s-wave component 
s = —0.25 are included to model the orthorhombicity of YBCO. Born impurity scattering with 

= 0.1 meV is also included. The bottom frame gives the critical frequencies Uc defining the van 
Hove singularities in the top frame. 

FIG. 7. Top frame: the impurity scattering rate ■~f{p, uo) in the unitary limit with = 0.1 meV 
as a function of the distance from the vortex core p, and frequency u. Here the magnetic energy 
v = 0.1, H II c-axis, and the vortex winding angle /? = 0°, i.e. the STM tip is moved along the 
antinodal direction. A system of tetragonal symmetry is considered. Bottom frame: the same 
as the top frame but for a system of orthorhombic symmetry having an effective mass anisotropy 
a = 0.4 and an s-wave component s = —0.25 to model YBCO. 
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